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ABSTRACT 



Aims. We reanalyse optical spectra of the z = 6.7 gamma-ray burst GRB 080913, adding hitherto unpublished spectra, in order to 
reassess the measurement of the neutral fraction of the IGM at high redshifts. 

Methods. In the data reduction we take particular care to minimise systematic errors in the sky subtraction, which are evident in 
the published spectrum, and compromise the analysis. The final combined spectrum has a higher S/N than the previously published 
spectrum by a factor of 1.3. 

Results. We find a single significant absorption line redward of the Lyor continuum break, which we identify with the SH+Sill 
AO. 126/jm blend, at z = 6.733. The sharp spectral break at Lyo- implies a comparatively low total column density of neutral hydrogen 
along the line of sight, log (iVm/cm ) < 20. We model the absorption with a host-galaxy DLA, surrounded by an ionised region of 
unknown size r, within the IGM of neutral fraction, x m . Despite knowing the source redshift, and the improved S/N of the spectrum, 
when fitting only over wavelengths redward of Lya, no useful constraints on *hi can be obtained. We consider the possibility of 
including the ionised region, blueward of Lya, in constraining the fit. For the optimistic assumption that the ionised region is trans- 
parent, Tap <K 1, we find that the region is of small size r < 2 proper Mpc, and we obtain an upper limit to the neutral fraction of the 
IGM at z = 6.7 of x m < 0.73 at a probability of 90%. 
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1. Introduction 

The hydrogen in the Universe became predominantly neutral 
at the epoch of recombination, z = 1070. The inter-galactic 
medium (IGM) out to z = 5, however, is highly ionised. Two 
observations made within the last decade have narrowed the 
window within which cosmic hydrogen reionisation took place. 
First, the electron scattering optical depth to the cosmic mi- 
crowave background measured by the Wilkinson Microwave 
Anisotropy Probe imp lies a redshift of reionisation of 1 1 .4 + 1.4 
1 dDunklev et al.ll2009l) . Second, observations of the Lya forest in 
the spectra of z ~ 6 quasars indic ate that z _ 5.8 marks the tail- 
end of the epoch of reionisation dFan et alj|2006l) . These two re- 
sults suggest that reionisation is an extended process. Therefore, 
in order to understand the chronology of reionisation in de- 
tail, there is considerable interest in detecting sources beyond 
z _ 6.4, the redshift of the most dis tant quasars so far discovered 
dFan et al.ll2003tlWillott et alJl200l . 

L ya emitting galaxies (LAEs) have be en detected out to 
z = 7 dKashikawa et alJ2006r.lbta et al.l2008h . Since the strength 
of the Lya emission line depends on the neutral fraction, xm, 
one approach advo cated to chart the progress of reionisation 
dMcOuinn et alJl2007l) is to measure the evolution of the proper- 
ties (e.g. abundance, clustering, etc) of LAEs. Howev er, the pro- 
cesse s involved are difficu lt to model accurately (e.g. Tasitsiomi 
2006; Diik straetal.ll2007h . and there is little consensus yet on 
the interpretation of the results from this approach. 

The confir mation of the co smological nature of gamma ray 
bursts (GRBs) dKulkarni et al.ll 19981) revealed their potential as 



probes of the high-redshift Universe. Bar kana & Loebl d2004l) 
emphasise two particular advantages of GRBs over quasars for 
studying reionisaton: that their afterglows will be visible to 
redshifts of z ~ 10; and that the size of the ionised region 
around the host galaxy, which complicates the interpretation of 
the spectrum, will be sma ll. Recent numerical simulati ons by 
iMcOuinn etafj d2008l) and lMesinger & Furlanettol d2008l) . how- 
ever, suggest that GRB host galaxies may be located in dense 
environments, and therefore the IGM immediately surrounding 
the host galaxy may be ionised by nearby quasars or local mas- 
sive star forming galaxies, and so the size of the ionised region 
needs to be included in the modelling. In any case a high signal- 
to-noise ratio (S/N) spectrum is needed to distinguish between 
the different signatures of a high-column de nsity of neutral hy- 
droge n in the GRB host galaxy (or nearby) ( R uiz-Velasco et all 
120071) . and of distributed neutral hydrogen in the IGM. 

Until the launch of the Swift satellit e dGehrels et al| | 2004 
the furthest known GRB was at z — 4.5 dAndersen et alJl200C 
comp ared to z = 5 .8 for the furthest known quasar at that 
time dFanetalJ koOO). Since then, three GRBs at z > 6 have 
been discovered using Swift : GRB 050904 at z = 6.3 in 
2005 dCusumano et al.ll2006l: iKawai et al.ll2006T) ; GRB 080913 
at z = 6.7 in 2008 dGreiner et al.ll2009h ; an d, in 2009, the re- 
markable source GRB 090423 at z = 8.2 dTanvir et al.ll2009L 
ISalvaterra et alj|2009l) . the most distant object yet found. 

iTotani et aH (|2006) present a detailed analysis of the opti- 
cal afterglow spectrum of GRB 050904, which was measured 
with the highest S/N of the three z > 6 GRBs. The spectrum 
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displays a red damping wing from Lya absorption that is pro- 
duced by some combination of a high-column density absorber 
near the GRB (hereafter referred to as a DLA, for damped Lya 
absorber), and a smoothly distributed component in the IGM. 
By fitting absorpti on models containing these two components, 
Tota ni et al.l (|2006) were able to place a limit on the neutral frac- 
tion of the IGM at z = 6.3 of x m < 0.17(0.60) at 68%(95%) 
confidence. The constraints are relatively weak despite the rea- 
sonably high S/N of the spectrum. The sources GRB 080913 
and GRB 090423 are potentially more interesting, because of 
their higher redshifts, but the published spectra have insufficient 
S/N to place any useful constraints on xhi when employing two- 
component (DLA+IGM) fits. 

Here we describe and analyse an improved spectrum of GRB 
080913, which includes unpublished spectroscopic data taken 
three nights after the published spectrum. In Section 2 we de- 
scribe the observations taken on each night, and the reduction 
techniques. We analyse the combined spectrum in Section 3, 
first searching for absorption lines in the new spectrum, in or- 
der to measure the redshift of the source, and then fitting a two- 
component DLA+IGM model to the observed continuum break. 
Finally, the results are summarised in Section 4. 

2. Observations and Data reduction 

iGreiner et al.1 (120091) (hereafter Paper I) provide a summary of all 
the photometric and spectroscopic observations of GRB 080913. 
Here we recap the details of the spectroscopic observations only. 
The object was first observed at the Very Large Telescope using 
the FOcal Reduction and low dispersion Spectrograph 2 instru- 
ment on the night beginning Septempber 13 2008. A 1 arcsec 
slit and the 600z grism were employed, providing wavelength 
coverage from 0.7470 to 1.0570jum. One 1800s and one 600s 
exposure were taken, but the second exposure is not useful as it 
was taken in twilight. The first spectrum is published in Paper I, 
and used for the analysis presented there. Due to poor weather, 
further observations were not possible until three nights later, on 
the night beginning September 16 2008, when 7 exposures of 
1800s each were secured. Because the source was then fainter, 
the S/N of the combined spectrum from the second night is lower 
than that of the single spectrum from the first night. 

We reduced the data in the following manner. After the stan- 
dard bias subtraction and flat-fielding steps, cosmic rays were re- 
moved from individual frames using t he Lap lacian Cosmic Ray 
Removal Algorithm of Ivan Dokkuml ((2001). The 7 frames ob- 
served on the second night were taken using 4 different slit posi- 
tions. We adapted methods from n ear-infrared spectroscopy de- 
scribed in IWeatherlev etal] d2005l) to minimise residuals from 
the subtraction of bright sky lines. After subtracting a functional 
fit up each column in every frame (first-order sky subtraction), 
all frames at other slit positions were averaged and subtracted 
(second-order sky subtraction). The 7 frames were then regis- 
tered and combined, weighting by the inverse variance. Since 
only one frame was taken on the first night, frames from the 
second night were used for the second-order sky subtraction. 
Compared to the standard method for optical long-slit spec- 
troscopy (i.e. first-order sky subtraction), the method followed 
here results in a slight increase in the random errors, but largely 
eliminates systematic errors. 

One-dimensional optimal spectral extraction was performed 
separately for each night, weighting by the profile of a nearby 
bright star. Wavelength calibration was achieved using observa- 
tions of a HeNeAr lamp. Corrections for telluric absorption and 
flux calibration were applied simultaneously using observations 



of the standard star LTT 7987. We scaled the spectrum from the 
second night to the spectrum from the first night, and combined 
the data, weighting by the inverse variance. The flux calibration 
agrees well with that in Paper I. 

3. Analysis of the Optical Afterglow Spectrum 

The final combined spectrum, binned to critical sampling, is 
shown in Fig.Q] The S/N of the new spectrum compared to the 
original spectrum shows an improvement by a factor of 1 .3 on 
average. The sharp break shows no clear evidence of a significant 
red damping wing, and is indicative of a relatively low column 
density of neutral gas, compared to GRB 050904. 

3.1. Search for absorption lines 

The Lya- continuum break near 0.94/mi corresponds to z ~ 6.73. 
In the hope of measuring an accurate redshift, we made a search 
for absorption lines in the afterglow spectrum. 

For all the fits presented in this paper we assumed that the 
spectral energy distribution of the GRB follows a power-law 
(F v oc v~ a ), where a = 1.12, which was derived from photo- 
metric measurements taken on the first night (Paper I). After 
subtracting the continuum fit, we tested for the presence of an 
absorption line centred at each pixel redward of the break, and 
then refined the measurement of the line centre of any detected 
lines, by refitting with the central wavelength as a free parameter. 
We assumed a Gaussian profile matched to the spectral resolu- 
tion (i.e. 8 unbinned pixels, FWHM = 0.00128/mi). The S/N at 
any pixel j is given by: 

(S/N)j = / 1 , (1) 

where u describes the line profile, centred on i = 0, f r is the 
residual flux, and <x is the flux error dBolton et al. 2003- 

Over the rest-frame wavelength range measured, one of the 
strongest absorption lines typically seen in afterglow spectra 
is the SH+Sill 0.1260/mi blend, which is expected to be sit- 
uated near 0.974/im. We found only one absorption line with 
significance greater than S/N = 2, detected at the wavelength 
0.9743/mi atS/N = 2.9. The line is marked in Fig.Q] The wave- 
length matches the expected location of the SH+Sill 0.1260^/m 
line, which is unlikely to be a coincidence. Therefore we infer 
that the line is real. The measured wavelength provides an ab- 
sorption redshift of z a b s = 6.733. 

3.2. Fitting the continuum break 

The continuum break is due to a combination of absorption by 
neutral gas in the GRB host galaxy and the IGM, therefore only a 
two-component (DLA+IGM) fit provides a meaningful physical 
model. The modelling is complicated by the unknown size and 
neutral fraction of the surrounding ionised region. For this rea- 
son it is useful first to fit single-component models to quantify 
the strength of the break. 

For each single-component model, the power-law SED is 
modified by absorptio n by either a DLA o r the IGM, computed 
using the equations in Totani et aQ (120061) . We limit the fits to 
data blueward of 1 .00/itri, and mask out the SH+Sill absorption 
line. For the DLA model we assume complete absorption blue- 
ward of the line centre, since this is achieved for a minimal value 
of xhi ~ 10~ 5 . The fits are determined using x 2 minimisation. 
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Fig. 1. Combined VLT/FORS2 spectrum of GRB 080913. The data have been binned by a factor 4, providing 2 binned pixels per 
resolution element. The l<x error spectrum is plotted in red and is offset by -2.0 x 10" 17 Wm~ 2 jum ~' for clarity. The position of the 
SH+Sill absorption line (as indicated), implies a source redshift of 6.733. The dotted line corresponds to redshifted Lya. 
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Fig. 2. Best-fit single-component DLA and IGM models, over- 
plotted on the GRB afterglow spectrum, which is rebinned by a 
factor of 4. 



The DLA model has three free parameters: the redshift, 
Zdla; the column density, Nm', and the continuum normalisa- 
tion at a wavelength of 0.96yum, cdla- We find best-fit values 
of log(A/ ffl /cm- 2 ) = 19.84, and zdla = 6.731. The redshift 
is in excellent agreement with the absorption-line redshift. The 
value of Nui is low in compa rison with the majo rity of GRB- 
DLA systems detected so far dFvnbo et al.l 12009]). which is in- 
teresti ng as this is in line with the prediction of Nagamine et alj 
(2008) that the typical column densities of GRB -DLA systems 
decrease towards higher redshifts. The IGM model is defined by 
four parameters: the neutral fraction xm, the continuum normal- 
isation at 0.96jum, cigm, and the upper and lower redshift lim- 
its over which x^i applies, zigm.u and Zigm.i- We fix the lower 
redshift to Zkjmj = 6.0. We find best-fit values of jchi = 0.21, 
and zigm.u = 6.737. The two models are plotted in Fig. [2] 
and are almost indistinguishable, implying that there will be a 
strong degeneracy between the parameters Nhi and xhi in a two- 
component fit. 

These two fits a re quite different from those obtained by 
iGreiner et alJ (2009), who found much stronger absorption. For 



the DLA model their 95.4% confidence range is 20.29 < 
log(AWcnr 2 ) < 21.41, and for the IGM model they found 
a best fit value x m = 1.00, with x m > 0.35 at 95.4% con- 
fidence. From an analysis of the two spectra, the differences 
may be explained by systematic sky-subtraction residuals in 
their spectrum, caused by the strong OH sky lines in the re- 
gion 0.93 - 0.96/mi (visible in the error spectrum, Fig. [TJ. We 
demonstrate this by a quantitative comparison of the two spectra 
and error arrays. We first binned each spectrum by a factor of 
four, summing the variance arrays appropriately. We then sub- 
tracted a median-filtered version of the data, and divided the 
result by the error array. In the absence of systematic errors, 
or significant features in the spectra, these S/N spectra should 
have mean zero and o"s/n ~ 1.0. We find this to be true for 
our spectrum , at all wavelengths, and also for the spectrum of 
IGreiner et al.l (|2009) in regions free of strong sky lines. However 
over the wavelength range of interest 0.93 - 0.96/im we measure 
cs/n = 1.7 in their spectrum, which is strong evidence of sys- 
tematic sky-subtraction residuals. 

Turning now to the two-component model, we firstly fix 
Zdla to the absorption redshift of z — 6.733, since the absorp- 
tion presumably arises in the dominant neutral-gas system in the 
host galaxy. However, we must consider the possibility that zdla 
and zigm.u differ, because the GRB host galaxy is surrounded 
by an ionised region. The two-component model therefore re- 
quires five parameters: xm, Zigm.u, Zdla, Nhi and the continuum 
normalisation at 0.96/vm, c. A further complication is the un- 
certain degree of absorp tion within the ionised bubble. The GP 
dGunn & Petersonl lT965) optical depth is given by 



r GP = 5.4 



1 +z 
7.7 



3/2 



+ A), 



(2) 



where A is the overdensity. If tcp is large, there is negligible 
transmission blueward of Lyes', so that only the wavelength range 
redward of Lye is useful. Under these conditions, at this low 
S/N, we find we can obtain no useful constraints on the various 
parameters of the two-component model - for the IGM com- 
ponent, which only contributes significantly to the wing of the 
absorption profile, a large ionised region and large x^i provides 
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Fig. 3. The posterior probability distribution of the neutral frac- 
tion when a joint DLA+IGM fit is employed. All other parame- 
ters have been marginalized. 

a similarly good fit as a small ionised region and small xm- On 
the other hand, if tqp within the ionised region is small, then, for 
a large ionised region, there will be measureable trans mission to 
the blue of the DLA line centre: see llliev et al.l (l2008l) for illus- 
trations of such cases from simulations. Under these conditions, 
with high S/N, this region could also be used in the modelling. 
For our spectrum, due to the low S/N, we consider only the sim- 
ple (optimistic) case that the ionised region is completely trans- 
parent, tqp <k 1. We assume a uniform prior on each parameter 
except for the column density. Since GRBs are most likely to 
occur in regions with significant neutral hydrogen, we apply a 
prior which scales linearly with column density. 

We find best-fit parameters of log (A^i/cirr 2 ) = 19.60, 
jchi = 0.06, Zigm.u = 6.733. To determine the uncertainties in 
each of these parameters we use a Markov Chain Monte Carlo 
(MCMC) algorithm and sample the posterior probability distri- 
bution. The result is presented in Fig.[3j which shows the poste- 
rior probability distribution of the IGM neutral fraction having 
marginalised over all other parameters. We find jchi < 0.73 with 
a probability of 90%. Our analysis also provides constraints on 
the size of the ionised region around the GRB host. In Fig. |U 
we plot the marginalized posterior probability contours in xm-r 
space produced from our fits. Here we see that an ionised re- 
gion which has a size smaller than ~ 2 proper Mpc is favoured. 
The regions in parameter space where r is negative correspond 
to situations where the DLA material is accelerated towards us 
due to either the GRB event or its progenitor. If we ignore these 
regions, we find that r < 1.3 proper Mpc with a probability of 
90%, which suggests that a large ionised region is not present 
around the host galaxy. This value is, however, still consistent 
with external source s of ionising flux being present. For compar- 
ison lHaimanl (|2002) found that the ionised regions around LAEs 
have a typical size of 0.8 proper Mpc. 

4. Conclusion 

New optical spectroscopic observations of GRB 080913 have 
been presented and analysed. The detection of SH+Sill absorp- 
tion (0.1260/mi) at 2.9cr provides a redshift of the DLA host 
galaxy of z = 6.733. Employing a joint DLA+IGM model to 
fit the observed continuum break, we find an upper limit to the 
neutral fraction of the IGM x m < 0.73 at a probability of 90%. 
However this result rests on the assumption that the ionised re- 
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Fig. 4. Posterior probability contours of the joint DLA+IGM fit 
in r-xm space. All other parameters have been marginalized. 
Regions enclosing 68.2%, 95.4% and 99.7% of the posterior 
probability are highlighted. 

gion surrounding the host galaxy is transparent to Lycr. Any 
analysis of a GRB spectrum needs to include the radius of the 
ionised region as a free parameter, and to consider the question 
of the neutral fraction within this zone. Furthermore the scat- 
ter in measurements of xhi between diff erent sources at simi- 
lar redshifts is predicted to b e substantial (McO uinn et alj [2008: 
iMesinger & Furlane tto 2008), and needs to be quantified. Higher 
S/N spectra of several sources at high redshift will be required 
to make significant progress in this field. 
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